The Statistics of Crumpled Paper 
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A statistical study of crumpled paper is allowed by a minimal ID model: a self-avoiding line bent 
at sharp angles - in which resides the elastic energy - put in a confining potential. Many independent 
equilibrium configurations are generated numerically and their properties are investigated. At small 
confinement, the distribution of segment lengths is log-normal in agreement with previous predictions 
and experiments. At high confinement, the system approaches a jammed state with a critical 
behavior, whereas the length distribution follows a Gamma law which parameter is predicted as a 
function of the number of layers in the system. 

PACS numbers: 68.55.-a,46.32.+x, 46.65.+g 



When a sheet of paper is crushed, the formation of 
a network of ridges is observed. The process is irre- 
versible as energy accumulates in small regions leading 
to localised plastic flow. Many fundamental questions 
arise. What is the resistance of the crumpled ball to me- 
chanical forces? What are the mechanisms of the cascade 
of energy to small scales? What is the probability distri- 
bution of lengths and energies? The two latter are similar 
to the questions central to hydrodynamic turbulence 0. 
Early experimental studies |2, LJ |J] of crumpled paper 
balls focused on fractal properties such as the scaling of 
the radius with the size of the flat sheet. The same prop- 
erty was used to characterise the phases of microscopic 
membranes - such as red blood cells or graphitic ox- 
ide - for which thermal fluctuations are relevant; these 
microscopic membranes raise a number of numerical and 
theoretical difficulties, particularly when self-avoidance 
is implemented 0. Crumpled paper was also viewed as 
a self-affine surface which roughness was predicted and 
measured [3, |£g ■ 

From elasticity theory 0, we know that thin elas- 
tic plates have two modes of deformation: bending - 
which involves curving the plate - and stretching - 
which changes the distances on the plate. Bending is 
much less expensive energetically than stretching so that 
pure bending deformations should be always preferred. 
However, this is not possible for a number of bound- 
ary conditions as shown in [inj : this leads to the for- 
mation of a near-singular network of lines (stretching 
ridges) and points (developable cones) where the ex- 
pensive stretc hing energy is localised. Even if isolated 
ridges [H 13 El or d-cones [IS H E El are rather 
well-characterised, the understanding of a full network is 
far from being complete. Most experimental and theo- 
retical studies tackled situations with a small number of 
singularities llil \vk Eil Eij. Eof or a highly symmetric 
network |2ll l22f . 

Our aim is to investigate numerically and theoreti- 
cally the statistics of crumpled paper. The strength 
of crumpled sheets was measured and found to involve 



logarithmic relaxation and a critical behavior close to a 
compact configuration [2^ |. Experiments also showed a 
broad distribution for the length of ridges, either indi- 
rectly - through noise emission [24 . 123 1 - or directly - 
through geometrical measurements This could be 

explained |27| in terms of the hierarchical random break- 
ing of ridges into smaller ones. However, full numerical 
approaches 0,113,12^1131 are impeded by the complexity 
of the problem and did not allow statistics on crumpled 
paper. Here we introduce a minimal ID model. We gen- 
erate numerically a large number of equilibrium configu- 
rations and study the resulting mechanical and geometri- 
cal properties. In particular, we show a transition of the 
length distribution from log-normal to Gamma and we 
predict the parameter of the Gamma distribution using 
arguments analogous to those used for mixing |3lT | and 
spray formation |32j . 

A 2D crumpled sheet can be considered as a self- 
avoiding polyhedron which energy is concentrated at the 
edges (stretching ridges |13|) and vertices (d-cones (lCf). 
The ID equivalent that we introduce is a self-avoiding 
closed line bent at sharp angles and which energy is con- 
centrated at the vertices (Fig. la). The important fea- 
ture here is that the position of the vertices along the 
line is free. This model can be viewed either as a rep- 
resentation of a cut through a crumpled ball of paper 
(Fig. lb) or as accounting for a crumpled sheet such 
that all ridges are parallel and have the same length 
L. The state of the system is defined by the coordi- 
nates (x-,y-) of the N vertices. The total length of the 
line is conserved and set to 1 by defining its energy as a 
function of the rescaled coordinates (xi,yi) = (x^y'^/l, 

I = EJ(4+i - <? + (24+1 " 24)T /2 - In order to keep 
the features of 2D sheets, the energy of each vertex is 
taken as the energy of a ridge (an edge) bent with the 
same angle <f>i - defined so that 4>i = when the sheet 
is flat - and of length L (the value of L is unimportant 
as it will be scaled out). The line is put in a confining 
quadratic potential of strength A. The self-avoidance is 
implemented through a hard-core interaction S^c with in- 
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FIG. 1: a A compact equilibrium state of the ID model. A 
self-avoiding line bent at sharp angles <j>i with N segments 
of free lengths 1% is put in a confining quadratic potential 
of strength A. Here the number of vertices is N = 80, the 
thickness h = 10 -3 (in units of total line length) and the 
confinement A = 10 6 . b Picture of a cut through a crumpled 
ball of paper, courtesy E. Couturier. 

finite cost (a large number in the simulation) between all 
pairs of segments, taking into account their thickness h. 
The total energy reads 

£ = E< tan 7/3 fa + A J r 2 (s)ds + £ hc . (1) 

r(s) is the (rescaled) position vector as a function of the 
curvilinear coordinate along the line. The unit of energy 
is the bending modulus of the sheet multiplied by the fac- 
tor (L/ft,) 1 / 3 that comes from the ridge energy [J^; the 
unit of length is the total length of the line. The number 
of vertices ./V and the thickness h - which appears only in 
the hard-core interaction - are left as parameters whereas 
the confinement A is the control parameter. This system 
has a large number of energy minima - which is favored 
by self-avoidance - and our aim is to investigate the sta- 
tistical properties of the corresponding equilibrium con- 
figurations. Note that we do not account for the plastic 
flow which reduces the stiffness of the ball and introduces 
a history dependance in the equilibrium states. 

In the numerical minimisation of the energy, Powell's 
algorithm [33| appeared as the most convenient to cope 
with the discontinuities of the hard-core interaction. This 



algorithm is built upon a set of M directions (M be- 
ing the number of degrees of freedom) which are used 
to determine the directions for the successive ID min- 
imisations. In order to obtain independent samples of 
the metastable states of the system, the set of directions 
is chosen at random. A better convergence was obtained 
by running Powell's algorithm 4 times and initialising the 
set of directions before each run. Care was taken in opti- 
mising the evaluation of the hard-core interaction which 
is expensive (typically 90% of CPU time) because it in- 
volves all pairs of segments. The results reported here 
correspond to 4 values of N x 6 values of the thickness 
h (in the range 10 -5 to 10~ 2 ) x 150 samples x 40 values 
of the confinement ~ 10 5 minimisations with 100 to 200 
degrees of freedom. 

We first measured the gyration radius R g = (r 2 ) 1 / 2 
averaged over all samples as a function of the confine- 
ment A (e.g. Fig. 2 for N — 80). For a given sample, the 
gyration radius depends on the choice of the sequence 
of applied confinements A = Aj., A2 . . . as in the experi- 
ments of . This can be ascribed to the large number of 
metastable states so that the equilibrium configuration is 
selected by the history of the system. The process of av- 
eraging over samples almost suppresses this dependence. 
At low lambda, the equilibrium state is a regular closed 
polyhedron so 2nR g ~ 1. When A reaches a threshold 
A c ~ 1, the polyhedron buckles and starts to fold. When 
A ~ 100, many contact points have formed. Eventually 
R g decreases slowly with A to a minimal value R g value 
such that the system is jammed; we roughly find 

J? g -4~A- a ) (2) 

with a = 0.3 ± 0.1. In this large A regime, we also found 
the total energy £ and the energy of the ridges £ r to scale 
as 

£~A /3 , £ r ~A ft , (3) 

with (3 = 1±0.05 and T = 0.45±0.05. The exponents can 
be understood as follows. The system reaches a jammed 
state of radius R g , so that the main contribution to its 

energy is the confining potential ~ \R3 and (3=1. The 
main contribution of ridges comes from U-turns, i.e. from 
the few ridges with angles </>i = 7r/2 — e, e <C 1 being a 
small angle. The radius of gyration approaches its lim- 
iting value as R g — R g ~ e 2 . Balancing the energy of 
ridges e~ 7 / 3 with the correction to the confining poten- 
tial - Ae 2 yields a = 6/13 ~ 0.46 and ft = 7/13 ~ 0.54. 
These exponents are slightly larger than in the numer- 
ics. The value of a differs from the experimental |23| 
and numerical [3(J value a = 0.54 for 2D sheets probably 
because of the different dimensionality. 

We examined the correlations of the normals along the 
line and we found them to be small except at high con- 
finement or for zero thickness as some pairs of segments in 
contact are exactly aligned. The probabilty distribution 
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FIG. 2: Gyration radius of the line R g as a function of the 
strength of the confining potential A. Continuous line: aver- 
age over 60 samples; dotted (resp. dashed) line: maximum 
(resp. minimum) value of P g among all samples. N — 80 and 
h = 10~ 3 . The differences between maximum and minimum 
are relatively large; nevertheless, the standard deviation is 
only of the order of 5 10 -4 . 

of angles is symmetric with respect to zero and almost 
flat - of course this probability vanishes for the maximal 
angles 4> = ±ir. In contrast, the distribution of distances 
between two consecutive vertices appeared to be broad. 
We first checked that a threshold on length (up to 10 h) 
or angle (up to 10°) - i.e. suppressing from the statistics 
small segments or vertices with a small angle - does not 
affect the results. 

At small confinement (l<A<10 4 orl> 2itR„ > 
0.2), we observe sequential buckling events such that 
a part of the line folds inward reducing the radius of 
the ball. As predicted in [27} and found experimentally 
in [2^, the hierarchical splitting of a line leads to a log- 
normal distribution (see Fig. 3a) with probability density 

Pln(z = t/(£» = ^exp (-^- + \) . (4) 

y/TTUX \ <7 Z 4/ 

This distribution is in agreement with our data (except 
for sizes comparable to the thickness h), with a width a in 
the range 1.0-1.8, comparable to the experimental values 
in the range 1.2-1.4 [26j. Note that these experiments 
were performed at low confinement - 2irR g > 0.6 in units 
of sheet length. 

In contrast, at high confinement (A > 10 4 or 2nR g < 
0.2), we found that a Gamma distribution with density 

P T ( X = £/{£)) = ^ expi-ax) (5) 
1 [a)x 

provides a better fit to the data than a log-normal dis- 
tribution which fails for the most likely sizes (Fig. 3b). 
This can be accounted for as follows. The buckling pro- 
cess is no longer hierarchical which at first sight should 
lead to an exponential distribution of lengths Pe(x = 
t/(£)) = exp(— x) as the result of the random splitting 
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FIG. 3: The distribution of segments lengths P (£/{£)) nor- 
malised by the mean length (N = 60 and h = 10" 3 ). a Low 
confinement A = 900. The line is a best fit to a log-normal 
distribution P LN (Eq.gJ with width a = 1.59 ± 0.03. b High 
confinement A = 10 6 . The two lines ares a best fits to a log- 
normal distribution Pln (Eq.[IJ and to a Gamma distribution 
Pr (Eq.EJ with parameter a = 1.64 ± 0.03. 

of a line |3^]- However, at large confinement, many self- 
contact points exist and segments in contact tend to cor- 
relate. Let us define the number of layers Al as the mean 
number of intersections of the system with a straight half- 
line starting from the center of the potential (see Fig. 1), 
averaged over all directions and all samples. Al varies 
from 1 at low confinement to about 5 at the largest con- 
finements (except in the case h = 0: A^l can be as large 
as 7.5). Note that the number of layers is always smaller 
than the simple estimate l/(27rPL g ) (in the range 1-8) 
which reflects the "roughness" of the ball. The system 
can be viewed as Al interacting layers. If they were 
independent, the distribution of lengths would be expo- 
nential Pe- Here the lengths of segments in contact tend 
to average. So the the length should be the mean of 
Al exponential random variables, i.e. a random variable 
with a Gamma distribution Pr of parameter a = Al. 
We measured the parameter from the fit to the numeri- 
cal distributions and we found that, as soon as Al > 2.5, 

a = a A L - b (6) 

with a = 0.95 ± 0.1 and b = 2.05 ± 0.2 (Fig. 4). This in 
agreement with our argument except that the effective 
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FIG. 4: The parameter a of the Gamma distribution Pr 
(Eq. |KJ as a function of the number of layers A^l (defined 
as the mean number of intersections of a half-infinite straight 
line with the system, averaged over all directions of the line 
and all samples). 



number of layer is smaller by 2 than Njj, which simply 
means that the system is slightly less correlated than if 
two segments in contact had exactly the same length. 
This robust behavior shows that the statistics is com- 
pletly determined by the average number of layers in the 
ball. 

To summarise, we have introduced a simple ID model 
which allows a comprehensive statistical study of crum- 
pled paper. At low confinement, a hierarchy of buckling 
events leads to a log-normal distribution of lengths in 
agreement with the predictions of and the exper- 
iments of |2fij . At high confinement, jamming is ap- 
proached with a critical behavior; self-avoidance leads to 
self-correlations and to a Gamma distribution of lengths 
solely determined by the number of layers in the system. 
This increase in correlations accounts for the Hurst expo- 
nent of compact paper balls || being larger than in the 
loose balls of j^. Obviously, our results call for experi- 
mental measurements of length distributions at high con- 
finement. Besides, our model can be refined in two ways. 
First, the number of vertices could be made free but one 
needs to find a criterion to add or suppress vertices. Sec- 
ond, it could be extended to 2D which is more realistic 
but yields a more difficult numerical task : if the sheet 
is assumed to be a polyhedron, developability becomes 
a local constraint at each vertex whereas self-avoidance 
becomes even more computationally expensive. 
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